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^ The  spatial  variation  of  many  measurable  quantities  such  as  pressure  or 
velocity  may  generally  be  synthesised  by  the  sum  of  a given  number  of  sinusoids. 
The  power  average  of  a small  number  of  random  samples  is  presented  as  a suitable 
estimate  of  the  mean  power  contained  in  such  a sum  of  sinusoids  of  random  argument 
and  equal  amplitude.  The  confidence  levels  which  apply  to  these  estimates  are 
presented,  and  the  technique  is  illustrated  by  its  application  to  the  measurement 
of  acoustic  radiated  noise  in  shallow  water  by  the  appropriate  use  of  vertical 
transducer  arrays. 
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INTRODUCTION 


In  many  practical  circumstances,  measurable  physical  quantities  such  a6 
pressure  or  velocity  display  significant  spatial  variations  which  may  be 
represented  by  the  sum  of  a given  number  of  sinusoids.  If  the  variation  is 


periodic  with  period  xq  the  spatial  variation  may  be  described  by  a discrete 


Fourier  series  of  fundamental  wavenumber  2lT/x  . Thus 

o 


p(x)  = a sintk  x + S ) 
m=1 


where 


= PJim/x 


If,  however,  p(x)  is  not  periodic,  the  fundamental  wavenumber  2U/x 


becomes  vanishingly  small  and  p(x)  may  only  be  described  by  a Fourier  integral 
which  represents  a continuous  spectrum  of  wavenumber  components.  Thus 


n 


w(k)e^kxdk 


In  this  case,  the  sinusoidal  components  are  harmonically  unrelated  and 


the  wavenumbers  k may  be  considered  to  represent  values  of  an  evenly  distributed 


random  variable. 


The  mean  power  contained  in  p(x)  may  be  defined  by  G^Cx)  where 


Gm(x)  = lim  i TP  p2(x  ) 

M N~*N  . n 

n=1 


where  the  x are  N random  locations  of  x. 
n 


Unfortunately,  in  cases  where  p(x)  defines  a spatial  variation,  it  is  usually 
impractical  to  collect  more  than  a few  samples  of  the  value  of  p(x)  since  each 
new  sample  represents  a new  transducer  location.  In  these  circumstances  an 
estimate  of  GH(x)  must  be  derived  from  a limited  number  of  samples.  If  thiB 

estimate  is  denoted  g^(x)  then 


6NM(x)  = 1 p2(xn) 


GM(x)  = N^^ogNM^X^ 


As  will  be  shown  later  in  this  note,  the  convergence  of  the  estimate  g.m(x) 

NM 


to  its  true  value  G^(x)  is  rapid  as  N increases,  enabling  an  accurate  estimate 


of  the  mean  power  to  be  made  using  as  few  as  ten  values  of  p(xn)  evaluated  at 


the  ten  random  locations  x 
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In  the  following  analysis,  all  of  the  contributing  sinusoids  are  assumed 
to  have  random  arguments  and  unit  amplitude.  This  latter  assumption  is  necessary 
to  render  the  problem  tractable.  In  these  circumstances,  it  will  be  shown  in 
Section  2 that  for  a given  number  of  samples  N,  the  confidence  which  may  be 
placed  in  the  estimate  g^Cx)  reduces  monotonically  as  M,  the  number  of  contri- 
buting sinusoids  increases.  For  the  special  case  M = 1 , it  will  also  be  shown 
in  Section  3 that  increasing  background  noise  levels  also  lead  to  significant 
monotonic  reductions  in  confidence  levels  for  signal  to  noise  ratios  of  less 
than  20dB.  For  an  infinite  number  of  contributing  sinusoids  of  equal  amplitude, 
confidence  levels  have  been  determined  by  application  of  the  chi-squared  distri- 
bution. In  Section  4,  the  applicability  of  the  confidence  levels  is  illustrated 
in  relation  to  the  enhancement  of  acoustic  noise  measurements  taken  in  shallow 
water,  by  the  appropriate  processing  of  outputs  from  transducers  which  form  a 
vertical  array. 


2.  CONFIDENCE  LIMITS  FOR  MEAN  POWER  DERIVED  FROM  RANDOM  SAMPLES  OF  THE  VECTOR 
SUM  OF  A VARIABLE  NUMBER  OF  PURE  SINE  WAVES  ‘ 

Consider  the  random  variable  p(x)  defined  by 

M 

p(x)  = 5^  sin(k  x) 
m=1  m 

where  the  wavenumbers  k are  taken  from  an  evenly  distributed  random  variable 

such  that  O^k  m 

m 

The  mean  square  value  of  p(x)  which  may  be  denoted  G^Cx)  is  given  by 
G„(x)  = 

due  to  the  orthogonality  of  the  trigonometric  functions. 

gNM^’  an  eetimate  the  mean  square  value  G^(x)  may  be  defined  as  the 

power  average  of  N samples  of  p(x)  evaluated  at  the  random  locations  x . 

n 


Thus 


snh(x)  - 

™ n=1 


V*>  is  related  to  the  estimate  g^Cx)  by  the  expression 

°m(x)  * 

The  accuracy  of  the  estimate  g^Cx)  may  be  defined  by  the  confidence  level 


that  any  estimate  gNM<x)  lies  within  the  bounds 

VX)  ^ gNM(x)«o<GM(x>,  (*>  1) 
o(. 

Lj^ot)  may  be  derived  from  the  expression 


J 


A 
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where  C^Cy)  i»  the  cumulative  probability  function  of  gNM(x)  such  that  C^Cy) 
is  the  probability  that  g^Cx^y. 

It  ia  also  convenient  to  define  Pj^x)  as  the  probability  density  function 
of  gJx)  since  Pm(x)  and  ^(y)  are  related  by  the  expression 


CNM(y)  =J°r  PNM(x)  *** 
o 

2.1  Derivation  of  the  confidence  levels 

As  shown  above,  the  confidence  levels  L-_,(«0  may  be  derived  directly 
from  the  appropriate  probability  density  function  P^(x).  Unfortunately 

the  derivation  of  P^(x)  cannot  be  carried  out  analytically  except  in  the 

special  case  in  which  both  M and  N are  1.  In  this  case  P..(x)  and  C (y) 
are  given  by  1 


V*> 


cnW 


ITx"  (1  - x)’ 
2 sin”1  (y^) 


For  M = 1 , values  of  P^(x),2^N<l6,  have  been  derived  from  P^(x) 

by  iterative  application  of  an  approximate  algorithm  which  was  designed  to 
calculate  joint  probability  density  functions.  A description  of  this 
algorithm  and  details  of  its  expected  accuracy  are  given  in  the  Appendix. 

For  values  of  M other  than  M = 1,  and  M =••,  values  of  Pj((M(x),1^N<l6, 

have  been  derived  solely  by  digital  simulation  on  a PDP  11/40  minicomputer. 
For  each  combination  of  N and  M such  that  1<N£l6,  and  M = 1,  2,  4,  8 or 
16,  100,000  values  of  the  estimate  g^M(x)  have  been  calculated  for  the 

assessment  of  confidence  levels.  The  derivation  of  PN<^(x),1#N^l6,  by 

digital  simulation  acts  as  a cross-check  on  the  accuracy  of  the  iterative 
calculation  algorithm  described  in  the  Appendix  and  on  the  convergence  of 
the  digital  simulation  process. 

For  M =«o,  the  amplitude  of  p(x)  must,  by  the  central  limit  theorem 
be  normally  distributed.  The  implication  of  this  is  that  gj^x)  must  be 

chi-squared  distributed  with  N degrees  of  freedom,  thus  enabling  values  of 
W°  to  be  derived  analytically  from  that  distribution. 

2.2  Numerical  Results 

The  normalised  values  of  PN<j(x)  and  values  of  CNl(x)  for  1^N^l6 

are  shown  in  Figures  1 and  2 respectively.  As  N increases  gNl(x)  tends 

to  a normally  distributed  variable  as  predicted  by  the  central  limit 
theorem.  In  Figures  3»  4 and  5*  values  of  L^(o0  are  shown  for  N and  M 
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such  that  1^N*l6  and  M = 1 , 2,  4,  8,  16  or»0.  The  values  ofet correspond 
to  tolerance  bounds  defined  by  G^(x)  ildB,  G^Cx)  -2dB  and  G^(x)  i3<*B. 

From  Figures  3,  4 and  5 it  can  be  seen  that  for  given  M and**,  ^(©0 

increases  monotonically  with  N.  In  addition,  for  given  and  N (N>l)  it 
can  be  seen  that  L^CoO  decreases  monotonically  as  M increases.  This 

second  conclusion  is  important  since  in  many  practical  situations  the  value 
of  M will  not  be  known  explicitly.  In  these  circumstances  however,  it  will 
be  possible  to  state  the  upper  and  lower  bounds  of  the  confidence  level 
that  may  be  applied  to  any  estimate  g^x),  (N>l)  since  for  given  N and  oC 

these  are  defined  respectively  by  L^Cot)  and 

In  particular,  if  G^(x)  -3dB  can  be  accepted  as  a tolerance  band,  then 

by  squaring  and  averaging  sixteen  samples  taken  at  random  from  any  known 
or  unknown  one-dimensional  distribution,  the  likelihood  that  the  average  is 
an  accurate  estimator  of  the  true  mean  square  of  the  distribution  will  lie 
between  94#  and  100#. 

Clearly,  one  conclusion  that  must  be  drawn  is  that  a far  better 
estimate  of  the  mean  power  contained  in  a spatial  measurement  environment 
may  be  obtained  by  power  averaging  data  derived  at  a relatively  small  number 
of  locations  than  may  be  derived  from  a single  measurement. 

. THE  EFFECTS  OF  GAUSSIAN  RANDOM  NOISE  ON  ESTIMATES  OF  MEAN  POWER  FOR  A SINGLE 


SINE  WAVE 

For  a single  pure  sine  wave  contaminated  by  Gaussian  random  noise  p(x)  may 
be  defined  by 

p(x)  = sin(kx)  + n(x) 

where  n(x)  is  the  instantaneous  level  of  noise  at  x. 

Since  M = 1,  G^Cx),  the  mean  power  contained  in  the  pure  signal  is  given  by 

G-,Cx)  = i 

Em.(x),  the  estimate  of  G (x)  is  defined  as 
N1  1 

8»i(x)  • £ibp2<v 

h=1 

where  the  x are  N random  values  of  x. 
n 


SN  as 


In  these  circumstances,  it  is  necessary  to  define  the  signal  to  noise  ratio 


SN  = 201og10  f rms  value  of  the  signal 
V rms  value  of  the  noise 


If  the^ noise  is  assumed  to  be  Gaussian  random,  then  SN  is  related  to  the 
variance  by 

S'2  = iio-SN/1° 

for  a unit  amplitude  sinusoidal  signal. 
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Thus,  for  example,  a signal  to  noise  ratio  of  OdB  may  be  simulated  by 
contaminating  the  pure  unit  amplitude  sinusoidal  signal  with  Gaussian  random 
noise  of  zero  mean  and  variance  of  0.5. 

The  effects  of  noise  on  values  of  L^C0*)  have  been  calculated  by  digital 

simulation  for  values  of  SN  ranging  from  OdB  to  40dB.  The  values  of  1^  («0 

already  shown  in  Figures  3,  4 and  5 apply  effectively  to  an  infinite  signal  to 
noise  ratio.  For  each  combination  of  N and  SN,  100,000  values  of  the  estimate 
g^(x)  have  been  calculated  in  assessing  the  noise  affected  confidence  levels. 

Figures  6,  7 and  8 show  the  confidence  levels  L^(®<)  that  the  estimate 
gN<J(x)  lies  within  the  tolerance  bands  G^xJlIdB,  G^Cx^dB  and  G^xJ^dB 
for  values  of  SN  equal  to  0,  10,  20,  3^»  40  andoo. 

As  might  be  expected,  the  introduction  of  increasing  levels  of  Gaussian 
random  noise  into  the  pure  sinusoidal  signal  causes  a monotonic  reduction  in 
confidence  levels,  although  this  degradation  is  insignificant  for  signal  to 
noise  ratios  in  excess  of  20dB. 


4. 


APPLICATION  TO  ACOUSTIC  NOISE  MEASUREMENT  IN  SHALLOW  WATER 


The  apparent  acoustic  source  radiation  measured  at  a point  in  a bounded 
medium  ir.  invariably  composed  of  the  phased  contributions  from  many  interfering 
propagation  paths  which  result  from  reflections  at  the  boundaries  of  the  medium. 
In  this  context,  shallow  water  may  be  defined  as  a bounded  medium  in  which  the 
effects  of  multipath  propagation  are  significant  at  any  required  measurement 
point . 

In  the  idealised  case  of  a fluid  laye’  bounded  by  a free  surface,  z = h 
and  a hard  surface  •£  = o,  the  far  field  acoustic  radiated  pressure  due  to  an 
omnidirectional  point  source  situated  at  (0,zq)  may  be  expressed  at  radian 

frequency  w by  a normal  mode  expansion  as  shown  in  Reference  1.  Thus, 
p(r,z,t)  = 2iwg  exp  (?)(W  £ cos  (b^zQ)cos  (bLz).  ^”^.exp(i  ^Lr)exp(iwt) 


where  b^  = (L-^)TT 

h 

and  h3L  = [(kh)2  - (L-J)^!2]^  L = 1,  2 .... 

(r,  6,  z)  are  cylindrical  coordinates  as  defined  in  Figure  10. 

It  is  convenient  to  define  a mode  number  M such  that  is  real  and 
- is  imaginary.  It  is  then  important  to  note  that  modes  defined  by  m 


(4.1) 


mode 


number  L<M  propagate  efficiently  to  the  far-field  whereas  modes  defined  by  L>M 

attenuate  exponentially  away  from  the  source.  Under  conditions  in  which  V 

1 

is  imaginary,  all  of  the  normal  modes  attenuate  exponentially  away  from  the 
source.  In  addition  for  L<M  the  power  propagated  by  the  efficient  normal  modes 
increases  monotonically  with  mode  number  L.  Thus,  the  far-field  acoustic  pressure 
may  be  expressed  by  the  sum  of  the  pressures  contributed  by  M normal  modes  each 
of  which  has  sinusoidal  depth  dependence. 
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In  general,  it  is  required  that  measurements  of  the  far-field  acoustic 

pressure  generated  in  shallow  water  by  a directional  source  should  yield  the 

mean  power  that  is  propagating  through  the  cylindrical  sector  defined  by 

r = r = constant  and  0 -S<e<9  + S where  & is  a function  of  the  source 
o o o 

directivity. 

In  sea  water,  suitable  measurements  may  be  made  by  a vertical  array  of 
randomly  spaced  transducers  which  are  located  at  the  positions  (r  ,9  ,z^ ),rQ 

and  9 being  constant, 
o 

For  an  omnidirectional  source,  situated  in  an  ideal  cylindrical  channel, 
the  far  field  acoustic  pressure  measured  at  radian  frequency  w by  the  ith 
element  of  such  an  array  is  defined  by  (4.1)  which  for  constant  r and  6 may 
be  simplified  to  give 


p(z.,t)  s aT(r  ) cos(bTz. )exp(iwt) 

1 T - " O w 1 


The  power  spectrum  of  the  time  variation  of  p(z^,t)  may  be  defined  by 

its  amplitude  s(z.)  which  is  given  by 

1 M 

s(z.  ) = wC  aT(r  ) cos(bTz.) 

i Ic1  L o Li 

If  the  true  mean  power  propagating  through  the  cylindrical  sector  r = r , 

9 -&<9<6  + £ is  given  by  Gu(z)  then  an  estimate  of  Gu(z)  may  be  defined  ° 
o o n M 

as  hMU(z)  where 

” N . 

hNM(z)  = I 21  8 (zn} 

N . 
n=1 

the  zn  being  randomly  selected  values  of  z. 

h|^(z)  is  related  to  G^(z)  by  the  expression 


Gm(z)  = lim  hNM(z) 

N-fc-oo 

Although  the  estimates  h^(z)  and  gjw(x)  bave  identical  forms,  the 

confidence  levels  L^W)  which  apply  to  the  estimates  gj^Cx)  nay  not  strictly 

be  applied  to  the  estimates  h^(z)  since  s(z)  is  synthesised  by  harmonically 

related  sinusoids  of  unequal  amplitude.  However,  it  will  be  demonstrated  that 
the  confidence  levels  Ljjfc)  and  L^aC)  may  be  considered  as  upper  and  lower 

bounds  for  the  confidence  levels  which  apply  to  values  of  h^(z). 

In  order  to  assess  the  applicability  of  to  estimates  h^Cz),  a 

further  digital  simulation  has  been  carried  out  to  derive  the  confidence  levels 
TNm(°0  which  may  be  applied  to  the  estimates  h^Cz)  for  the  particular  measure- 
ment geometry  in  which  a horizontal  range  of  100m  exists  between  the  measuring 
array  and  a source  situated  at  mid-depth  in  water  6lm  deep.  For  this  horizontal 
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range,  the  exponentially  decaying  modes  will  still  contribute  significantly 
to  the  measured  radiation  in  the  frequency  range  0-30Hz  which  is  equivalent 
to  values  of  M42.  The  application  of  the  confidence  levels  T^(«0  and  T^  {■<) 

to  measurements  using  this  geometry  will  thus  be  erroneous  in  practice. 
However,  the  quantity  T^(®0  - will  give  a good  measure  of  the  effects 

on  confidence  levels  of  nor  -uniform  normal  mode  amplitudes  and  of  harmonic 
relationships  between  the  individual  propagating  normal  modes. 

In  the  derivation  of  T^(o<)  each  value  of  M effectively  defines  a 

frequency  range  in  which  exactly  M propagating  normal  modes  exist.  For  each 
M,  this  frequency  range  has  been  split  into  ten  thousand  bands,  and  at  each 
band  centre  frequency  f an  estimate  hj^(z)  has  been  formed  by  averaging  the 

squared  pressures  calculated  from  (4.1)  at  N randomly  generated  depths.  The 
confidence  levels  T.^(«<)  have  been  derived  by  comparing  each  value  of  h^(z) 

with  the  respective  true  mean  square  pressure  propagating  at  the  frequency 
f . For  each  combination  of  N and  M,  only  10,000  values  of  h._.(z)  have  been 

computed  due  to  the  time  consuming  nature  of  the  calculation  described  by 
(4.1). 


Values  of  T^oO  are  presented  implicitly  in  Figure  9 which  gives  values 
of  - Lj^**).  Since  TN^(o<)  and  L^(o0  should  be  identical  for  given  N 

and  °<,  the  slight  deviations  of  T^(o<)  • L^(®<)  from  zero  give  a measure  of 
the  inaccuracies  introduced  by  the  digital  simulation  processes. 

From  Figure  9,  it  is  apparent  for  N>3  that  is  strictly 

positive  for  M>1 . For  other  combinations  of  N and  M,  the  quantity  " 

V00  is  generally  positive  except  in  a few  isolated  cases  in  which  the 

negative  excursions  are  within  the  bounds  that  might  be  expected  to  result 
from  inaccuracies  inherent  to  the  digital  simulation  procedures.  It  nay  also 
be  deduced  from  Figure  9 that  for  N>3,  is  always  bounded  by  the  confidence 

levels  (°<)  and  L^c*).  Thus 

^ tnm(°°  * W°°  (N>5’  1^M^0> 

This  is  an  extremely  important  conclusion,  since  it  implies  that  L^O*) 

and  Lj^0*)  act  as  bounds  for  the  confidence  levels  that  apply  to  estimates  of 

mean  radiated  power  that  have  been  derived  in  shallow  wacer  by  power  averaging 
the  frequency  domain  outputs  of  randomly  spaced  elements  that  form  a vertical 
array. 

For  an  array  of  12  randomly  spaced  elements,  the  level  of  confidence  that 
the  power  average  of  the  individual  array  element  outputs  is  within  ±3dB  of 
the  true  mean  radiated  power  lies  between  89«5%  and  99. 5%.  The  equivalent 
confidence  bounds  for  an  array  of  16  elements  are  94%  and  100%. 


4.1  Experimental  results  of  noise  propagation  trials,  undertaken  in 
shallow  water 


In  February  1975 » AMTE  carried  out  a series  of  noise  propagation 
trials  in  shallow  water  in  which  the  noise  radiated  by  a hydrosounder 
at  a horizontal  range  of  100m  was  measured  by  a vertical  five  element 
hydrophone  array.  The  hydrosounder,  which  is  practi  ally  omnidirectional 
in  the  frequency  ranpe  0 - 1kHz,  was  situated  at  a depth  of  J>0w  in  water 
of  depth  80m.  A single  hydrophone  was  situated  10m  directly  above  the 
hydrosounder  in  order  to  monitor  the  source  output.  This  monitor  hydro- 
phone may  be  considered  as  being  in  the  acoustic  far-field  of  the  source 
at  frequencies  above  300Hz.  The  layout  of  the  experimental  facility  is 
shown  in  Figure  10. 

Experimental  data  gained  from  these  trials  are  shown  in  Figures  11, 

12  and  l4.  All  of  these  data  have  been  analysed  using  Fast  Fourier 
Transform  techniques  as  implemented  in  Reference  3* 

Figure  11  shows  the  frequency  power  spectrum  of  the  monitor  hydro- 
phone response  to  a hydrosounder  output  of  the  form  of  shaped  band 
limited  random  noise  in  the  frequency  band  250Hz  - 1kHz.  The  nominal 
bandwidth  is  4Hz.  The  response  is  not  completely  smooth  due  to  the 
interfering  effects  of  the  first  surface  reflection.  From  Figure  12, 
which  shows  the  transfer  functions  between  each  array  element  and  the 
monitor  hydrophone,  it  is  clear  that  the  measured  radiated  pressure  may 
be  described  accurately  by  the  contributions  from  only  two  propagation 
paths.  These  are  the  direct  path  and  the  path  which  embodies  a single 
surface  reflection.  On  this  assumption  of  radiation  into  a half  space, 
calculated  theoretical  predictions  of  the  data  shown  in  Figure  12  ore 
shown  for  comparison  in  Figure  13.  The  agreement  is  remarkably  good  for 
the  array  elements  situated  at  depths  of  24,  JO  and  38  metres,  but  for 
the  elements  at  47m  and  56m,  the  effects  of  reflections  from  the  sea  bed 
cause  some  deterioration  in  the  detail  agreement. 

As  was  shown  earlier  in  this  section,  a good  estimate  of  the  mean 
power  radiated  by  the  hydrosounder  at  a range  of  100m  may  be  obtained 
by  averaging  the  power  spectre  cJnr.I.vcad'  from  the  individual  array  elements. 
However,  if  the  narrow  band  radiated  pressure  is  to  be  related  to  a known 
excitation, 3dB  must  be  added  to  the  spatial  average  of  the  transfer  function 
between  the  array  elements  and  the  excitation,  which  is  defined  by  the 
monitor  hydrophone.  This  is  because  the  narrow  band  source  spectrum  will 
always  describe  the  peak  square  value  of  the  excitation  in  each  narrow 
frequency  band.  Since  the  spatially  averaged  narrow  band  transfer  function 
should  relate  the  mean  square  radiated  pressure  to  the  mean  square  level 
of  the  excitation,  3dB  should  be  added  to  the  averaged  transfer  function 
in  order  to  convert  the  peak  square  level  of  the  excitation  to  its  more 
appropriate  mean  square  level. 

Figure  14  shows  the  result  of  adding  3^B  to  the  power  average  of  the 
experimental  transfer  functions  shown  in  Figure  12.  The  data  in  Figure 
14(a)  have  been  likewise  derived  from  the  theoretical  data  presented  in 
Figure  13.  The  major  conclusion  that  may  be  drawn  from  these  Figures  is 
that  the  power  averaging  of  even  a small  number  of  array  element  outputs 
leads  to  significant  improvements  in  the  accuracy  of  radiated  noise 
transfer  functions. 
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In  the  case  of  acoustic  radiation  by  a point  source  into  a half 
space,  p(z),  the  variation  of  pressure  with  depth  (z)  at  a fixed 
horizontal  range  r , is  not  explicitly  sinusoidal.  However  for  practical 
measurements,  p(z)°may  either  be  expanded  into  a discrete  harmonic  series 
or  into  a continuous  wavenumber  spectrum.  In  either  case,  for  a given 
number  of  array  elements,  N,  and  for  a given  tolerance  bandog  , a lower 
bound  given  by  Lj^oO  may  be  stated  for  the  confidence  levels  that  are 

applicable  to  a transfer  function  that  has  been  derived  by  power 
averaging  the  N frequency  power  spectra  of  the  individual  array  elements. 
Values  of  may  be  derived  from  Figures  3»  4 or  5« 

Figure  5 shows  that  the  use  of  5 array  elements  for  estimating  the 
mean  propagating  power  implies  that  the  likelihood  of  the  estimate  being 
within  ±3dB  of  the  true  mean  lies  between  70#  and  94#.  This  is  borne  out 
by  Figure  l4  which  shows  a scatter  of  ±4dB. 

Figure  15  shows  third  octave  transfer  functions  between  the  individ- 
ual array  elements  and  the  monitor  hydrophone.  Figure  16  shows  the  results 
of  adding  3dB  to  the  power  average  of  these  transfer  functions.  These 
Figures  also  include  theoretical  losses  for  point  source  radiation  into 
an  infinite  medium  and  into  a half  space.  From  these  Figures,  it  is 
apparent  that  third  octave  smoothing  is  useful  in  removing  the  effects  of 
propagation  whenever  the  source  spectrum  is  relatively  flat.  It  must  be 
stressed  however  that  for  other  source  characteristics,  third  octave 
averaging  is  not  suitable  for  removing  the  effects  of  multi-path  propaga- 
tion. 

From  Figure  16,  it  is  also  clear  that  the  addition  of  3dB  to  the 
averaged  transfer  function  provides  a suitable  estimate  of  the  mean  power 
propagating  away  from  the  source  at  a range  of  100m  since  the  propagation 
loss  between  the  source  and  the  array  is,  as  might  be  expected  of 
radiation  into  a half  space,  about  JdB  less  than  the  equivalent  free- 
field  propagation  loss. 

5.  CONCLUSIONS 

Confidence  levels  L^(o<)  applicable  to  g^(x)  have  been  derived  where 
®NM^  an  estimate  the  mean  power  contained  in  p(x),  the  sum  of  M sinusoids 
of  equal  amplitude  and  random  argument.  gNM(x)  and  p(x)  are  defined  by 

M 

p(x)  = S sin0rn 
m=1 

gNM(x)  = j£p2(xn) 
n=1 

where  0 and  x are  random  variables  with  even  distributions, 
m n 

L^CoO  is  defined  as  the  level  of  confidence  that  any  estimate  gNM(*)  lies 
within  lOlog^^dB  of  the  true  mean  square  value  of  p(x)  which  is  given  by  ^M. 
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It  has  also  been  shown  that  (®()  and  I^JoQ  act  as  bounds  for  the 

confidence  levels  which  apply  whenever  the  sinusoids  that  synthesise  p(x) 
have  unequal  amplitudes  and  are  harmonically  related.  This  enables  the 
to  be  applied  to  measurements  of  acoustic  radiated  noise  taken  in  bounded” 
media.  In  this  context,  it  has  been  shown,  both  theoretically  and  experiment- 
ally, that  significant  improvements  in  estimates  of  the  mean  radiated  power 
propagating  in  shallow  water  may  be  achieved  by  power  averaging  the  frequency 
spectra  obtained  from  the  individual  randomly  spaced  elements  of  a vertical 
array.  It  is  concluded  that  the  use  of  12  elements  in  this  way  leads  to  an 
estimate  that  is  between  90%  and  99%  certain  of  being  within  ±3dB  of  the  true 
mean  radiated  power.  Moreover,  the  confidence  levels  represent  lower 

bounds  of  the  confidence  levels  that  may  be  applied  to  estimates  where 

measurements  are  taken  of  any  type  of  signal  variation. 


A W Walker 
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FIG.  5 CONFIDENCE  LEVELS  Lnm(«)  THAT  THE  ESTIMATE  g^M  (X)  IS  WITHIN  ±3dB 
OF  THE  TRUE  SIGNAL  MEAN  SOUARF  VALUE  WHEN  N SAMPLES 
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FIG.  7 CONFIDENCE  LEVELS  THAT  THE  ESTIMATE  gNI  (X)  IS  WITHIN  i2dB 
OF  THE  TRUE  SIGNAL  MEAN  SQUARE  VALUE  WHEN  N SAMPLES 
ARE  TAKEN  RANDOMLY  FROM  A SINGLE  SINE  WAVE  CONTAMINATED 
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FIG.  9 DIFFERENCE  IN  CONFIDENCE  LEVELS  ( TNM  {«■ ) - L NM  («0 ) 

TOR  1*N*16,  M = 1, 2,  4,  8 AND  16  AND  FOR  VALUES  OF  «c 
EQUIVALENT  TO  TOLERANCE  BANDS  OF  MEAN  SQUARE  *1dB.  i2dB.  s 3dB 
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HYDROPHONE  RADIATED  NOISE  RESPONSE  TO  HYDROSOUNDER  OUTPUT 
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PROPAGATION  LOSS  BETWEEN  MONITOR  HYDROPHONE  AND  INDIVIDUAL  ARRAY  ELEMENTS  (THIRD  OCTAVE  SPECTRUM  LEVELS) 
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FIG.  15  EXPERIMENTAL  THIRD  OCTAVE  PROPAGATION  LOSSES  BETWEEN 
MONITOR  HYDROPHONE  AND  INDIVIDUAL  ARRAY  ELEMENTS  COMPARED 
WITH  EQUIVALENT  THEORETICAL  FREE -FIELD  AND  HALF-SPACE  LOSSES 
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FIG.  16  EXPERIMENTAL  THIRO  OCTAVE  PROPAGATION  LOSS  BETWEEN 
MONITOR  HYDROPHONE  AND  POWER  AVERAGED  ARRAY  OUTPUTS  *3dB  COMPARED 
WITH  EQUIVALENT  THEORETICAL  FREE  -FIELD  AND  HALF- SPACE  LOSSES 


APPENDIX 


CALCULATION  ALGORITHM  FOR  DETERMINATION  OF  JOINT 
PROBABILITY  DENSITY  FUNCTIONS 


The  calculation  algorithm  described  in  this  Appendix  has  been  developed 
specifically  to  allow  PN<)(x)  to  be  derived  by  iterative  means  from  P^Cx)  as 
the  sole  input. 

Strictly  P^(x)  should  be  derived  from  P^Cx)  alone  by  evaluation  of 
the  expression 

Pfj^  (y)  = •••• 

v s 

V x2+  ....  +xN  = y 

Numerical  evaluation  of  P^^y)  by  this  method  using  a digital  computer 

is  however  excessively  time  consuming  and  the  alternative  approximate  algorithm, 
which  is  the  subject  of  this  Appendix,  has  been  used  for  the  iterative 
generation  of  P^Cy)  from  P^(x)  and  x). 

It  is  convenient  to  descibe  this  algorithm  in  relation  to  the  derivation 
of  P^(y)  which  is  the  probability  density  function  of  the  average  of  the  squares 

of  N samples  taken  at  random  from  the  same  pure  sine  wave.  Any  particular  value 
of  the  sum  of  the  squares  of  N random  samples  taken  from  this  sine  wave  will  be 
denoted  by  ^(N)  where 

vf»(N)  =2~8in2& 

n=1 

0 is  taken  from  0 which  is  an  evenly  distributed  random  number  in  the  range 

°<0n4  Tr/2. 

The  range  of  ^(N)  is 

0^(N)<N 

In  order  to  combine  P(;g_1)1(y)  with  P^Cy)  to  produce  pNl(*)  it  is 

necessary  to  split  the  ranges  of  *4^(1 ) into  NB  bands  of  equal 

width.  The  use  of  500  such  bands  gives  very  good  convergence  in  this  case. 

PNl(y) , which  is  a continuous  function,  may  thus  be  represented  discretely 

by  NB  terms  of  the  form  j^.l-^m^NB,  where  ^P^  is  defined  as  the  probability 

that  ^(N)  will  fall  in  band  number  m.  The  range  of  values  of  ^(N)  covered  by 
jjP^  is  defined  by 

(m-l)Nh^4#(N)^  mNh  (h  = VNB) 


p11(x1)P11(x2)....P11(xN)dx1dx2....dxN 
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and  P^Cy)  is  related  to  by  the  expression 

- V(",-i)NhKNh 

To  create  from  P^Cy)  811(1  P(N-^)1^y^>  tlie  following  procedure  has 
been  adopted. 

The  range  of  values  of  tp(N)  covered  by  ^P^  is  known  to  be 
(r-l)Nh^v|>(NX  rNh 

For  each  n,  1 ^ n ^ NB,  the  range  of  values  of  ^(N-l)  covered  by 

„ .P  is  known  to  be 
N-1  n 

(n-l)(N-l)h*vJ/(N-l)^  n(N-l)h 

and  for  each  m,  l^m^NB,  the  range  of  values  of  vj/(i)  covered  by  ^P^  is 
known  to  be 

(m-l)h  mh 

If  the  two  probabilities  ^ Pr  and  ^P^  are  multiplied  together,  the  range 
of  values  of  *^(N)  covered  by  the  product  is  given  by 

(n-l)(N-l)h  + (m-l)h'£'fKN)<f  n(N-l)h  + mh 

If  this  range  overlaps  the  range  of  „P  by  an  amount  X , then  the 

n r inn 

contribution  to  ..P  from  the  combination  of  „ P and  ,P  is  defined  as 
IT  r N-1  n 1m 


X 

mn 

rNh-(r-1 )Nh 


*N-1Pn'lPm 


If  the  value  of  X is  negative, 

X is  set  to  zero.  10,1 
mn 

This  process  is  repeated  for  all 


that  is  the  two  ranges  do  not  overlap, 

possible  values  of  n,  m and  r yielding 

X « „P  ‘.P 
nm  • N-1  n 1 m 

Nh 


Successive,  iterative  application  of  this  algorithm  thus  enables  PN  (y) 
to  be  calculated  for  any  value  of  N using  P^y)  as  the  only  input. 

A.1  Checks  on  calculation  algorithm  accuracy. 

In  order  to  establish  the  validity  of  the  calculation  algorithm  both 
theoretical  and  practical  checks  have  been  carried  out. 
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The  theoretical  check  was  performed  by  using  the  calculation 
algorithm  to  generate  iteratively  the  chi-squared  distribution  for  N 
degress  of  freedom  using  as  the  only  input,  the  chi -squared 

distribution  for  one  degree  of  freedom  in  the  range  0<X2<20.  For 
each  value  of  N,  the  calculated  distribution  has  been  compared  with 
the  theoretical  chi-squared  distribution  for  N degrees  of  freedom 
given  by 


Figures  A1  and  A2  show  typical  comparisons  between  theoretical  and 
calculated  values  for  N = k,  16  respectively.  The  agreement  between  the 
theoretical  and  calculated  data  is  remarkably  good  since  this  check 
applies  stringent  conditions  in  that  the  range  of  values  for  chi-squared 
is  strictly  infinite  whereas  the  calculation  algorithm  can  only  accept 
data  in  a finite  range  of  values. 

As  a practical  check  of  the  calculation  algorithm,  values  of  P^(x) 

have  also  been  derived  by  digital  simulation.  Figures  A3  and  kb  show 
typical  comparisons  between  data  derived  by  the  two  methods  for  N = 3,  16. 
The  minimal  lack  of  agreement  is  due  to  slight  uneveness  in  the  random 
number  generator  used  for  the  digital  Simulation. 


PN(X2)  . 


H<!>  2iK 


FIG.  A1  COMPARISON  OF  THEORETICAL  PROBABILITY  DENSITY  FUNCTION  OF  CHI-SQUARED  DISTRIBUTION 
FOR  4 DEGREES  OF  FREEDOM  WITH  THE  PROBABILITY  DENSITY  FUNCTION  CALCULATED  FOR  THE 
SAME  DISTRIBUTION  USING  THE  ITERATIVE  ALGORITHM  DESCRIBED  IN  APPENDIX 


